201903.00083v1 


chinaXiv 


兰 志 勇 ” 男 1980 年 
生 ， 博 士 ， 主 要 研究 方向 
为 大 功率 永 磁 同步 电机 设 
计 及 优化 、 en 
设计 及 优化 、 精 密 伺服 驱 
动 器 的 研究 与 as 


魏 雪 环 ” 女 1990 年 
Et ， 硕 士 研 究 生 ， 研 
究 方 向 为 永 磁 同 步 电 机 及 
异步 电机 的 温度 场 分 析 及 
优化 等 。 


[ny 


站 


ChinaXiv 合 作 期 刊 
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摘要 : 永 磁 同步 电机 与 其 他 类 型 电机 相 比 性 能 指标 尤为 突出 ， 被 广泛 应 用 于 各 个 
领域 。 随 着 电机 制造 业 的 发 展 ， 电 机 单机 容量 增加 ， 使 电机 内 部 温 升 过 高 ， 影 响 电 机 
的 性 能 与 可 靠 运行 。 本 文 在 考虑 高 速 电机 永 磁 体 涡流 损耗 对 电机 温度 场 的 影响 的 情况 
下 ， 用 集 总 参数 热 网 络 法 (LPTN) 对 永 磁 同步 电机 进行 温度 场 分 析 ， 建 立 38 市 点 的 
热 网 络 模型 ， 真 实 反映 了 电机 各 部 件 的 温 升 。 最 后 ， 通 过 与 有 限 元 法 (FEM) 分 析 结 
有 果 、 电 机 温 升 试验 结果 对 比 ， 证 明 该 热 网 络 计算 模型 的 正确 性 。 

关键 词 : 永 磁 同 步 电 机 ”温度 场 ” 集 总 参数 热 网 络 法 ”有 限 元 法 
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Thermal Analysis of PMSM Based on Lumped Parameter 
Thermal Network Method 


Lan Zhiyone Wei Xuehuan LiHuru Liao Keliang Chen Linhong 
(Xiang Tan University Xiangtan 411105 China) 


Abstract: The performance targets of PMSM are prominent, compared to other type 
motors, PMSM are widely used in various fields. With the development of the motor 
manufacturing industry, the motor capacity increases, resulting in increasing temperature 
of inside the motor, the high temperature severely affected performance and reliable 
operation of motor. In this paper, in consideration of the influence of the permanent 
magnet eddy current loss of high speed motor on the temperature field, thermal analysis 
using LPTN in PMSM, the thermal network model of 38 nodes is established. truly reflect 
the temperature rise of each part of the motor. Finally, the correctness of the model is 
verified by comparison with the results of finite element method (FEM) and temperature 
rise test. 
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1 引言 


近年 来 ， 由 于 永 磁 同 步 电机 具有 体积 小 、 容 量 
大 且 效 率 高 等 优势 ， 在 各 个 领域 得 到 了 广泛 的 应 用 。 
随 着 电机 功率 密度 的 提高 和 电机 材料 的 充分 利用 ， 
电机 经 常 处 于 电 负 荷 和 磁 负 荷 的 极限 状态 ， 导 致电 
机 的 损耗 变 大 ，3 引 起 电机 的 温 升 过 高 。 对 于 永 磁 电 
机 而 言 ， 过 高 的 温度 容易 使 永 磁体 发 生 不 可 逆 退 磁 ， 
影响 电机 的 性 能 和 可 靠 运行 "。 因 此 ， 对 电机 温度 
场 的 分 析 尤为 重要 。 

目前 ， 常 用 的 电机 温 升 计算 方法 有 有 限 元 法 
(FEM)、 集 总 参数 热 网 络 法 (LPTN) 和 流体 力学 
(CFD)。FEM 和 CFD 都 是 数值 计算 方法 ， 可 以 通 
过 计算 机 分 析 得 到 电机 整体 温度 的 分 布 云图”。 但 
是 ， 对 于 当前 的 计算 机 设备 用 这 两 种 方法 计算 ， 耗 
时 也 很 长 。LPTN 法 则 通过 将 电机 相似 部 位 以 节点 
的 形式 划分 为 几 个 部 分 ， 借 助 电路 的 思想 ， 通 过 热 
阻 及 温度 计算 公式 ， 得 到 电机 内 部 各 节点 的 温度 。 
由 于 热 阻 的 计算 方法 相似 ，LPTN 对 于 不 同 的 电 
机 模型 也 能 做 相应 的 简化 和 拓展 。 相 较 于 FEM 和 
CFD ， 该 方法 市 省 了 很 多 时 间 ， 计 算 时 间 上 具有 较 
大 优势 。 

LPTN 在 感应 电机 和 同步 电机 的 温度 场 分 析 及 
其 优化 过 程 中 都 有 应 用 ,文献 [各 用 LPTN 对 电机 
进行 了 温度 场 分析 ， 本 文 在 其 基础 上 添加 了 6 个 端 
盖 节 点 、 两 个 轴承 节点 和 3 个 转轴 节点 ， 使 电机 的 
网 格 结构 更 加 规则 ， 更 符合 电机 传 热 的 实际 情况 。 
并 在 考虑 永 磁 体 测 流 损耗 的 基础 上 用 LPTN 对 电机 
进行 温度 场 分 析 与 计算 。 最 后 ， 与 FEM 分 析 结 果 
及 电机 温 升 试验 结果 对 比 ， 验 证 该 方法 的 有 效 性 。 


2 集 总 参数 热 网 络 法 (LPTN) 


LPTN 是 将 电机 的 主要 部 件 按照 其 各 自 结构 及 
发 热 散热 情况 划分 成 多 个 温度 单元 ， 并 通过 热 导 将 
各 个 单元 联系 起 来 ， 使 电机 内 部 温度 场 离散 为 等 效 
热 网 络 ， 编 写 Matlab 计算 程序 ， 对 电机 各 部 分 温度 
进行 求解 ， 得 到 电机 的 平均 温 升 分 布 情况 。 

2.1 分 析 过 程 

用 LPTN 对 电机 进行 温度 场 分 析 计 算 时 ， 需 要 
对 电机 的 计算 模型 进行 合理 的 假设 ， 如 下 : 

(1) 电机 的 温度 分 布 周 向 对 称 ， 同 时 认为 电机 
在 圆周 方向 的 冷却 条 件 相 同 。 

(2) 电机 内 部 两 侧 的 空 腔 中 ， 各 点 的 空气 温度 
相同 ， 计 算 时 可 以 用 一 个 点 计算 。 
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(3) 忽略 定子 槽 部 绕组 的 趋 肤 效应 。 

(4) 两 个 而 点 间 的 互 导 相等 ， 并 与 温 升 无 关 。 
(5) 忽略 电机 的 辐射 传 热 过 程 。 

LPTN 的 温度 场 分 析 过 程 如 图 1 所 示 。 
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图 1 LPTN 热 分 析 过 程 
Fig.l Thermal analysis process of LPTN 


2.2 热 网 络 模型 

LPTN 的 热 网 络 模 型 是 二 维 模 型 ,不仅 可 以 浪 
虐 轴 向 热 传 递 ， 也 可 以 考虑 径 向 热 传 递 。 

本 文 以 十 二 极 三 相 永 磁 同 步 发 电机 为 分 析 实 
例 ， 其 具体 参数 为 : 额定 功率 320kW， 额定 转速 
6 000rvmin， 定 子 外 径 0.4m， 定 子 内 径 0.28Sm;， 转 
子 外 径 0.282m， 转 子 内 径 0.Im， 电 机 模型 如 图 2 所 
示 。 电 机 轴 向 简化 模型 及 各 部 件 间 的 热 阻 分 布 如 图 3 
所 示 。 


图 2 电机 模型 
Fig.2 Motor model 


由 图 3 的 电机 各 部 件 轴 向 热 阻 分 布 关系 可 知 ， 
电机 的 每 个 相 邻 部 件 都 通过 热 阻 有 相应 的 联系 。 热 
阻 是 传导 热 阻 和 对 流 热 阻 ， 为 
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式 中 ，! 为 两 节点 间距 离 ，4 为 接触 面 面 积 ，4. 为 对 
流 截面 积 。 


3 ”电机 轴 向 热 阻 分 布 


Fig.3 Distribution of axial thermal resistance 


由 上 述 分 析 ， 对 样机 轴 向 分 布 进行 网 格 划 分 ， 
在 各 区 域 中 心 设立 节点 ， 使 网 格 剖 分 保持 整齐 ， 保 
证 计算 精度 ， 要 充分 考虑 电机 的 结构 、 材 料 和 工 亏 
等 因素 。 沿 轴 向 把 定 转子 划分 成 均等 三 份 ， 样 机 的 
热 网 络 模 型 节点 分 布 如 图 4 所 示 。 其 中 ，a-e,， 5 个 
市 点 为 电机 机 过 表面 市 点 、1-5 为 机 过 市 点 、6-8 为 
定子 轿 部 闻 点 、9-13 为 绕组 节点 、14-16 位 定子 此 
部 节点 、17-19 为 永 磁体 节点 、20-22 为 转子 节点 、 
23-24 为 机 腔 空 气节 点 、25-30 为 转轴 节点 〈 节 点 30 
为 轴承 节点 )、31-33 为 前 端 盖 节 点 、34-36 为 后 端 
盖 节 点 、37-38 为 轴承 节点 。 


全 一 OO-O-O- -9 一 G: 
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一 热传导 “和 旬 对 流散 热 
人 有 源 节 点 @ 无 源 节点 


4 热 网 络 模型 节点 分 布 


Fig.4 Node distribution ofthermal network model 
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2.3 热源 分 析 


2.3.1 热源 计算 

热源 与 电机 温度 场 的 求解 密 不 可 分 ， 它 与 电机 
的 各 部 分 损耗 有 关 。 永 磁 同 步 电机 的 损耗 ， 主 要 包 
括 定子 铁 耗 、 绕 组 铜 耗 和 机 械 损 耗 ， 这 些 损耗 都 可 
根据 经 验 公 式 进行 计算 。 一 般 情况 下 ， 大 多 数 
的 电机 温度 场 分 析 忽 略 了 永 磁 体 涡 流 损耗 对 温度 场 
分 布 的 影响 ， 如 文献 [7] 虽然 对 永 磁体 涡流 损耗 进 
行 了 分 析 ， 但 并 没有 分 析 永 磁体 涡流 损耗 对 电机 温 
度 场 的 影响 。 然 而 ， 对 于 对 温度 较为 敏感 的 永 磁体 
材料 而 言 ， 永 磁体 涡流 损耗 所 导致 的 温 升 是 不 能 忽 
略 的 。 

根据 二 维 有 限 元 涡流 损耗 分 析 方 法 ， 在 时 域内 
磁场 方程 为 


+VH. (3) 


Es 
ot 


1 
V-(VA4,)=J -0o 
u 


式 中 , / 为 相对 磁 导 率 ; 4z 为 磁 位 矢量 ;7 为 电流 

密度 ; o 为 材料 导电 率 ; 5 为 电势 标量 ，H 为 永 磁 

体 矫 奖 力 。 
绕组 中 的 电流 密度 

d4 


-一 -VE 


(4) 


J=-0 


当 磁 场 发 生变 化 时 ， 会 在 导电 材料 内 部 产生 感 
应 电流 ， 该 电流 为 涡流 。 滴 流产 生 的 损耗 即 涡 流 损 
耗 为 


_ rfl7z| 
a 2 
式 中 ,为 轴 方 向 的 电流 密度 分 量 。 
2.3.2 热源 分 布 

(1) 绕组 铜 耗 分 布 。 铜 耗 分 为 模 部 和 端 部 损耗 ， 
其 取 值 与 槽 部 和 端 部 的 绕组 长 度 成 正比 。 


上 (S$) 


覃 部 损耗 的 表达 式 为 
了 
人 (6) 
端 部 损耗 的 表达 式 为 
Ls 
fowa = Fo, 7 (7) 


N 


式 中 , 工 为 定子 铁心 长 ;Lz 为 绕组 半 古 长 ; Ls 为 绕 
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组 端 部 伸 出 端 。 

(2) 定子 铁 耗 分 布 。 本 文通 过 Ansoft Maxwell 
软件 对 电机 此 部 和 力 部 磁 密 进行 分 析 ， 如 图 5 所 示 。 
根据 丽 部 和 力 部 的 质量 比值 ， 可 以 计算 得 到 样机 此 
部 与 拒 部 的 铁 耗 ， 见 表 1。 


1 000 
800 
一 
E 600 
mm 400 
200 
0 
0 20 40 60 80 100 120 
D/mm 
(a) 斩 部 磁 密 
2.0 
1.6 
i 
SN 0.8 
0.4 
0 
0 10 20 30 40 50 60 70 80 90 
D/mm 
(b) 齿 部 磁 密 
5 定子 碰 窗 
Fig.3 Stator flux density 
表 1 定子 铁 耗 分 布 
Tab.1 Iron loss distribution of stator 
定子 罗 部 损耗 /W 定子 齿 部 损耗 /W 总 损耗 
1 369 1 241 2610 


2.4 热 阻 计算 

用 LPTN 对 电机 进行 温度 场 分 析 计 算 过 程 中 ， 
最 关键 的 因素 在 于 各 个 市 点 间 热 阻 的 计算 ， 它 直接 
关系 到 整个 计算 结果 的 准确 性 。 各 个 节点 是 否 有 热 
阻 连 接 ， 主 要 由 电机 各 部 位 的 散热 路 径 决 定 。 电 机 
散热 部 分 主要 包括 机 充 散 热 、 定 子 散 热 和 转子 散 
热 ， 文 献 [9] 中 介绍 了 简单 的 对 流传 热 和 热传导 计 
算 公 式 。 

一 般 情况 的 LPTN 热 模 型 忽略 了 电机 的 端 盖 市 
点 ， 假 设 电机 轴承 产生 的 热量 为 零 ， 即 认为 其 损耗 
为 零 。 但 对 于 高 速 电机 而 言 ， 轴 承 损耗 较 大 ， 忽 略 
轴承 损耗 将 会 增 大 电机 温度 场 计算 的 误差 。 因 此 ， 
在 以 往 LPTN 热 模型 的 基础 上 ,添加 6 个 端 盖 市 点 
和 2 个 轴承 节点 ， 更 能 反映 电机 传 热 的 实际 情况 。 
2.4.1 端 盖 散热 

(1) 市 点 31 的 热 阻 计算 。 节 点 31 与 端 盖 节 点 
32、 和 轴承 市 点 37 和 机 外 空气 进行 热传导 。 与 节点 
32、 市 点 37 为 简单 的 热传导 。 
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二 (23) 
2A.53137 24as93la7 
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式 中 ，B, 为 轴承 厚度 ; 1。 为 轴承 导热 系数 ，Zae 为 
端 盖 厚 度 ，Xas 为 端 盖 导热 系数 。 

(2) 而 点 32 的 热 阻 计算 。 端 盖 市 点 32 与 机 腔 
内 空气 节点 23、 端 盖 节 点 31、 市 点 33 存在 热 交 
换 ， 还 与 电机 外 部 空气 存在 对 流散 热 。 此 处 只 介绍 
与 节点 23 的 热 阻 。 


在 -7-27 To 
4A, S13 283223 


Ry 到 


(24) 


式 中 ，Ls 为 机 过 长 度 。 
2.4.2 轴承 散热 

以 布点 37 为 例 ， 市 点 37 与 转轴 节点 25、 机 
过 闻 点 31 存在 热 交 换 。 

节点 37 与 转轴 节点 25 的 热 阻 为 


D, ,D.-D, 
2 人 ii93725 24.S3755 


3725 一 (25) 
2.5 热平衡 方程 

根据 上 述 分 析 ， 通 过 电机 内 部 各 节点 间 的 热 传 
递 关 系 ， 可 以 得 到 电机 所 有 节点 的 热平衡 方程 组 ， 
其 矩阵 形式 为 


GT =W (26) 

式 中 ,，G 为 38 阶 热 导 和 矩阵; 7 为 38x1 温 度 列 矩 
阵 ， 球 为 38 x1 损耗 列 矩 阵 。 
2.6 LPTN 计算 程序 

在 得 到 电机 各 节点 间 的 热平衡 方程 后 ， 需 用 
Matlab 编写 LPTN 计算 程序 ， 主 要 包括 原始 数据 输 
入 、 热 阻 计 算 、 温 度 和 矩阵 计算 和 数据 输出 ， 从 而 得 
到 各 节点 温度 。 对 于 不 同 的 电机 模型 ， 可 以 快速 地 
修改 部 分 程序 语言 ， 得 到 电机 的 温 升 分 布 ， 体 现 了 
LPTN 的 简易 性 和 拓展 性 。 


3 ”样机 温 升 计算 结果 对 比 


文献 [10] 用 FEM 对 样机 温度 场 进行 计算 ， 并 
分 析 了 永 磁 体 涡 流 损耗 对 电机 温度 场 分 布 的 影响 。 
可 得 到 FEM 和 LPTN 对 样机 温度 场 分 析 的 结果 ， 
见 表 2。 

由 表 2 可 明显 看 出 ， 在 永 磁体 温度 的 计算 上 ， 
两 种 方法 的 计算 结果 误差 最 大 ， 主 要 在 于 气 隙 处 理 
方法 的 不 同 。FEM 将 气 阶 等 效 为 导热 系数 一 定 的 固 
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表 2 计算 结果 表 3 结果 对 比 
Tab.2 Calculation result Tab.3 Comparison results 

平均 温 升 LPTN FEM 误差 (%) 平均 温 升 LPTN FEM 测试 结果 
定子 思 /CC 86.924 85.126 2.1 绕组 /C 90.7 89.1 91.5 
定子 齿 /C 89.289 87.617 J:9 机 元/C 39.6 38.8 40.2 

绕组 /C 90.653 89.051 1.8 
永 磁体 /'C 50.147 48.208 4.0 S 结论 
转子 粥 /'C 48.273 46.836 3:] 

轴承 /C 二 i 本 文 用 LPTN 对 电机 进行 了 温度 场 分 析 ， 并 与 
商 善 /C ee ee jy FEM 温度 场 分 析 结 果 进 行 了 对 比 ，FEM 可 以 得 到 
机 志 /C 让 5 5 电机 的 整体 分 布 去 图 ，LPTN 可 以 得 到 电机 各 部 件 


体 材 料 ， 并 通过 等 效 计 算 公 式 计 算 气 隙 的 等 效 导 热 
系数 ， 而 LPTN 则 通过 经 验 计 算 公 式 计算 出 流 一 固 
间 的 对 流 热 阻 ， 然 后 进行 温 升 计算 。 因 此， 由 于 对 
气 隙 处 理 方式 的 不 同 导 致 永 磁体 的 散热 情况 不 同 ， 
也 直接 影响 了 转子 轿 部 的 温度 分 布 情况 。 但 在 误差 
允许 的 范围 内 ， 有 限 元 法 和 集 总 参数 热 网 络 法 的 计 
算 结果 一 致 。 
4 样机 温 升 试验 

320kW 样机 测试 实物 如 图 6 所 示 。 为 了 验证 
FEM 法 和 LPTN 法 计算 结果 的 正确 性 ， 需 对 样机 
进行 温 升 测试 。 将 5 个 热 敏 电阻 预先 坦 入 5 个 定子 
槽 内 ， 对 样机 的 绕组 温 升 进行 测试 ， 用 红外 线 温度 
传感器 对 样机 机 过 进行 测 温 ， 可 以 得 到 机 过 的 平均 
温 升 。 


图 6 样机 测试 实物 图 
Fig.6 Prototype test 


由 表 3 中 的 定子 绕组 和 机 这 平均 温 升 可 知 ， 在 
一 定 的 误差 允许 范围 内 ， 用 FEM 和 LPTN 对 电机 
进行 温度 场 分 析 都 是 合理 的 。 


的 平均 温 升 ， 二 者 结果 一 臻 。 通 过 与 电机 温 升 试验 
结果 对 比 可 知 ，FEM 和 LPTN 对 电机 进行 温度 场 
分 析 都 是 合理 的 ， 都 可 在 工程 计算 中 应 用 。 与 FEM 
相 比 ，LPTN 节省 了 很 多 时 间 ， 对 于 不 同 的 电机 模 
型 也 能 做 相应 的 简化 和 拓展 ， 在 计算 时 间 和 灵活 应 
用 上 具有 较 大 优势 ， 为 下 一 步 基 于 电机 过 热点 的 优 
化 设计 做 好 了 准备 。 
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